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It is demonstrated that a two-component Bose-Einstein condensate (BEG) with all-repulsive inter¬ 
atomic interactions loaded into a radially symmetric harmonic trap supports robust non-coaxial 
vortices with approximately orthogonal vortex lines in each of the components. These cross vortices 
are excited from the linear modes by a sudden switch-on of the nonlinearity (via Eeshbach resonance) 
and are characterized by persistent dynamical regimes of precession with nutation, resembling the 
motion of a rigid body. The obtained dynamics can be understood qualitatively on the basis of a 
simple mechanical model. 
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Quantized vortices are fundamental topological objects 
responsible for a variety of interesting phenomena ob¬ 
servable in such areas as superconductivity [T], superflu¬ 
idity [2], and Bose-Einstein condensates (BECs) [3l |4]. 
In a three-dimensional (3D) one-component BEG, where 
a vortex is characterized by its charge, as well as by the 
geometry and dynamics of the vortex line, vortices re¬ 
ceived considerable attention during the last two decades. 
Multiple studies were devoted to vortex shapes and their 
stability [SHin]- More recent investigations of BECs con¬ 
centrated on the dynamics of vortex lines, in particular 
Kelvin waves m , vortex clusters [12] , and dipoles [13] , as 
well as on the existence of more complex topological ob¬ 
jects like hybrid vortex solitons la and twisted toroidal 
vortex solitons m in inhomogeneous media. 

Even richer topological structures appear in multicom¬ 
ponent condensates m- These are, in particular, mix¬ 
tures of hyperfine states, like the ones explored since the 
very first experiments reported observation of vortices in 
BECs [T7142Q] . Such BECs have a spinor nature and are 
characterized by macroscopic order parameters with two 
or more components. Nowadays, the complex dynam¬ 
ics of the phase separation [21] as well as sophisticated 
topological states like vortex solitons [22] , skyrmions [23] , 
vortex lattices [24], and Dirac magnetic monopoles [25] 
in spinor BECs have been described theoretically and ob¬ 
served experimentally. 

Here we present non-coaxial vortices existing in two- 
component BECs, which have a significantly simpler ge¬ 
ometry than the topological defects mentioned above. 
The objects we report here consist of two vortex lines, 
each one belonging to one of the components of a bi¬ 
nary mixture (we term them cross vortices). Families 
of cross vortices solutions bifurcating from two orthogo¬ 
nal vortices existing in the linear decoupled condensates 
and differ in many aspects form the most of topological 
structures studied, so far. In particular, their genera¬ 


tion does not require trap rotation or linear coupling be¬ 
tween the components. Such cross vortices are observed 
in miscible configurations [26] , where each component has 
nonzero density both inside and outside the core region 
of the other component and is characterized by a hole 
where the vortex lines intersect. The reported objects 
can be excited by switching-on the nonlinearity and sup¬ 
port a plethora of persistent dynamical regimes when vor¬ 
tices interact being oriented by arbitrary relative angles. 
Among all such regimes, we concentrate here on persis¬ 
tent precession and nutation which closely resemble the 
dynamics of a heavy top [27] . 

We consider a binary mixture of BECs of the atomic 
components having equal atomic masses m, and loaded 
in radially symmetric parabolic traps, identical for both 
components and characterized by the linear oscillator fre¬ 
quency ujQ. We also assume that the intra-species scat¬ 
tering lengths ai ^2 are equal, i.e., ai = a 2 = a, and the 
inter-species interactions are characterized by the s-wave 
scattering length ai 2 . Further, we adopt dimensionless 
units where time and coordinates are measured in 2/ujq 
and ao = ^Jh/muj{) units, respectively. The mixture is 
described by the two-component macroscopic wave func¬ 
tion (spinor) 4^ = (' 0 i,' 02 )^ which solves coupled dimen¬ 
sionless Gross-Pitaevskii (GP) equations 

^ i’l (la) 

^ + ^^^2 + (|^/’2p + ?/)2 (lb) 

Here a = 012 !cl G [0,1] characterizes the relative strength 
of the inter- and intra-species interactions, i.e., we con¬ 
sider positive scattering lengths when the homogeneous 
mixture is miscible (a stratified phase is absent [26]). In 
the chosen normalization the total number of atoms is 
given hy M = Nao/{S7ra) with N = J 

Equations 0 conserve the number of atoms in each 
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component, i.e. A^i ,2 5 the total energy E = 
where Ej = j dx is the en¬ 

ergy of the j-th component and Eint = <^ J 
is the energy of inter-species interactions, as well as the 
total angular momentum L = Li -1- L2 where Lj = 
-if Ip* (r)(r X V)'ipj(r)dr. Meantime, neither the ener¬ 
gies Ej nor the angular momenta Lj are conserved as 
function of time. In particular one obtains 

^ =-^ = IV’ 2 p(r X V)|i/;ipcir, (2) 

i.e., the change of the individual components of the an¬ 
gular momenta occur because of the inter-species inter¬ 
actions. 

Usually, the generation of single vortices is achieved 
either by rotating traps [4] or by dynamical m and 
topological pniiiH] phase imprinting. Nonlinear inter¬ 
actions of our model Q allow us to adopt the following 
methodology. Initially two (linear) vortices are created in 
the absence of two-body interactions and thus indepen¬ 
dently in each of the components. Then the nonlinearity 
is switched on (by means of the Feshbach resonance). 
Even though this represents a rather strong initial ’’per¬ 
turbation” of the nonlinear vortex states, we demonstrate 
below that the objects created in such a way are robust 
cross vortices, which are accompanied by rotational dy¬ 
namics. We notice that growing a vortex state departing 
from the linear limit was previously discussed in the lit¬ 
erature [6], however, in the context of a one component 
condensate and for a gradual (rather than instant) in¬ 
crease of the inter-atomic interactions. 

Let us now turn to the stationary solutions of 0 with 
the components having equal chemical potentials /i, for 
which 2 pj = e“*^^i4j(r) with Uj{T) solving the stationary 
GP equations {j = 1,2) 

jjiUj = -V‘^Uj +r‘^Uj -h {\uj\‘^ a\us-jf) Uj. (3) 

We start with the linear eigenvalue problem = 

(—-hr^)^ (hereafter we use tilde to specify the eigen¬ 
values and eigenfunctions of the linear limit). The com¬ 
ponents are now decoupled and the eigenstates are given 
by the well known eigenmodes of the linear harmonic os¬ 
cillator. Without loss of generality, we fix the z-axis along 
the vortex line of the simplest linear vortex of the first 
component, i.e., 'Ui(r) = ^ which 

corresponds to the chemical potential /i = 5 and has Ni 
atoms (change of the direction of the vortex lines along 
axes will affect only the direction of the rotations dis¬ 
cussed below). The vortex line in the second component 
can be taken to be rotated by an angle P with respect 
to the second chosen axis. We chose it to be the ^-axis 
which means that U 2 {v') = + iy')e~^ 

where 

( cos/d 0 — sin P \ 

0 10. 
sin P 0 cos P j 




FIG. 1: (Color online) Isosurfaces of the components of the 
cross vortex corresponding to \uip — 1 (left upper panel) and 
\u 2 p — 1 (right upper panel) for fi — If) and dependence N 
vs fi for several values of the interspecies interactions. 


Gonsidering now the nonlinear family bifurcating from 
the spinor (fii,fi 2 )^ one can compute the interaction en¬ 
ergy which is deposited into the system when weakly non¬ 
linearity is switched-on: Eint{P) ~ (cos^ yd-h 3) . 

Thus EintiP) achieves its minimum at P = 7 r/ 2 , i.e. 
when the two vortex lines are mutually orthogonal, as 
illustrated in the upper panel of Fig. A family of 
nonlinear cross vortices [characterized by the dependence 
A/'(/i)] bifurcates from the state (fii,'U 2 )^ corresponding 
to p = 7 r/ 2 . The isotropy of the system with respect to 
the y axis suggests that the entire family, i.e., the non¬ 
linear cross vortices are also characterized by orthogonal 
vortex lines. Such families are shown in the lower panel 
of Fig. for several values of the interspecies interaction 
a. It is evident from this figure that all these families 
grow from a single solution corresponding to /i = /i = 5. 
Since there is no qualitative difference between the cases 
with different a we consider a = 0.5 hereafter. 

Turning now to the dynamical problem we simulate nu¬ 
merically the evolution of a cross vortex after an instan¬ 
taneous switch-on of the nonlinearity. We evolve the so¬ 
lution of 0 starting from initial conditions = 0 ) = 
2 . 12 (x-h ' 02 (t = 0 ) = 2.12{z-\-iy)e~^‘^ (the 

number of atoms in each component is Ni = N 2 ^ 25). 
Snapshots of the evolution of the binary condensate are 
shown in Fig.|^(a). We observe that both vortex lines ex¬ 
perience persistent rotations returning to the initial state 
after the full period T 248. In Fig.j^b) we present the 
evolution of the projections of Li^ 2 - The dynamics of 
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the angular momentum components shows two types of 
oscillations. The oscillations with a longer period and 
larger amplitude correspond to rotation of the individual 
momenta Li ^2 with respect to the conserved total angu¬ 
lar momentum L, i.e., the precession of the cross vor¬ 
tex. The precession is accompanied by oscillations with 
a smaller period and a smaller amplitude. The latter 
describe oscillations of Li^ 2 , whose amplitude is charac¬ 
terized by the angle 0 measuring deviations of Li ^2 from 
their initial orientations Lj(t = 0 ) = IjQj [see Fig. [^d) 
below]. The dynamics of this second type corresponds 
to nutation. Fig.j^c) shows a vectorial representation of 
the evolution of the angular momenta of each component 
Li,2- 

Our simulations demonstrate the robustness of the 
nonlinear cross vortices, even though they are excited 
by initial conditions significantly deviating from station¬ 
ary nonlinear vortices. The dynamics is accompanied 
by modulations of the vortex lines (and consequently of 
vortex shapes). Numerical studies performed for differ¬ 
ent N result in qualitatively similar regimes, however, 
the dynamical characteristics are significantly affected by 
the two-body interactions. This is illustrated in Fig. 
where we show the precession [Fig. [^a)] and nutation 
[Fig.j^b)] periods and the nutation amplitude [Fig.j^c)] 
vs the number of particles. For small N the precession 
and nutation are slow and they accelerate as N increases, 
while the nutation amplitude decreases. 

The obtained dynamics is not obvious for at least two 
reasons. First, earlier observed precession of vortices 
other types was explained by instabilities miH]. In our 
simulations, however, we did not observe a transient pe¬ 
riod describing development of instabilities, but the ro¬ 
tational dynamics started immediately after switching- 
on the nonlinearity. On the other hand, Eq. ([^ indi¬ 
cates that a ’’force” necessary to induce precession and 
nutation should be related to anisotropy of the atomic 
distributions in each of the vortices and cannot be in¬ 
duced at t = 0 by the imposed initial conditions [the 
right hand side of 0 is zero for any vortex of the form 
'ijjj = {x ^ iy)(j)j{r\^ where = [x^y)]. To un¬ 

derstand the observed evolution we use a simple mechan¬ 
ical analogy shown in Fig. (d). The angular momenta 
Loi and L 02 of the linear modes used for the initial con¬ 
ditions, can be viewed as an initial perturbation of the 
angular momenta of the nonlinear cross-vortex with Li 
and L 2 . Since now Li • L 2 7 ^ 0 , i.e., /3(t = 0 ) < 7 r/ 2 , the 
dynamics shows a precession whose period is faster for 
a smaller /3(t = 0) [see Fig. [^(b)]. Concomitantly, due 
to the inter-atomic interactions there appears a ’’force” 
tempting to restore /d = 7r/2 [i.e. to minimize Eint{P)]- 
This force is responsible for the nutation. We also ver¬ 
ify that ^(Li • L 2 ) 7 ^ 0 (notice that the components of 
Li and L 2 oscillate out-of-phase) which means that the 
modula Li ^2 are not constant. We did not found sta¬ 
tionary vortices with (3 < 7 r/ 2 , but only precessing ones 



FIG. 2: (Color online) (a) Snapshots of the cross vortex dy¬ 
namics with N=25 at the instants of time indicated in the 
panels. The isosurfaces correspond to |'0i,2|^ = 0.1. Blue and 
red colors represent the and '02 components, respectively, 
(b) Solid (dashed) lines show the dynamics of the x (blue), y 
(red), and z (green) components of the angular momentum Li 
(L2). The moduli of Li and L2 are shown by the overlapping 
purple lines and the conserved modulus of the total angular 
momentum L is shown by the solid black line, (c) The evolu¬ 
tion of Li (left panel) and L 2 (right panel). The single green 
lines represent the invariant L. Circles with arrows indicate 
axes positions at t = 0 as in upper left snapshot of panel (a). 

[Fig. El (b)]. This means that P < it j 2 defines ’’initial” 
angular velocity of the precession (/3 = it/2 corresponds 
to zero velocity). This difference in the initial angular ve¬ 
locities is clearly visible when comparing panels showing 
the dynamics in Fig. for different initial angles. 

Both precession and nutation disappear in the linear 
limit {N 0): the periods become infinite. Thus 1/N 
can be viewed as a small parameter for the estimate of 
the periods. Figs.j^a) and (b) show the respective fitting 
curves. The dependence of the angular nutation ampli¬ 
tude is approximated by an almost linear dependence on 
N^ see Fig. [^c). 
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FIG. 3: (Color online) (a) Precession periods obtained nu¬ 
merically (circles) and interpolated by = 142.6 + 

2824.228/iV — 2.375/iV^ (dashed line), (b) Nutation periods 
obtained numerically (circles) and interpolated by = 

3.5 + 160/A^ (dashed line). (c) Nutation amplitude ob¬ 
tained numerically (circles) and interpolated by = 

9.5 + 0.056iV - 5.2 • - 2.3 • (degrees, dashed 

line), (b) Schematic diagram of the geometry of the angular 
momenta. 



Although, the presented ’’mechanical” picture is over¬ 
simplified, it results in reasonable estimates for the dy¬ 
namical parameters. If {Lj) is the average value of 
|Lj|, the average angle of precession can be estimated 
as 0' = arccos(2(I/i)/I/) and the average angle of nuta¬ 
tion as ^ = 7r/4 —respectively. Thus the ratio between 
the nutation and precession periods can be estimated as 
Tnut/Tpr = sin(^')/sin(^). Using these simple estimates 
for the case of = 25 shown in Fig. we obtain 0' « 
0.75196, 0 ^ 0.03344 and hence TnutiTpr ~ 20.42. The 
result of our simple considerations is not very far from 
the numerically observed value of (Tnut/Tpr )~ 25.05 
but clearly deviates from it. This difference can be ex¬ 
plained by the strength of the inter-atomic interaction 
which results in a significantly perturbed nonlinear dy¬ 
namics. Reducing the number of particles to = 10 (not 
shown here), i.e., reducing the influence of the nonlinear¬ 
ity, one obtains much ’’cleaner” precession and nutation. 
For this case, our estimate yields 0' « 0.7537, 0 ^ 0.0317 
and hence Tnut/Tpr ~ 21.6 while numerically we find 
(Tnut/Tpr )~ 21.9 which is a remarkable accuracy. 

Remarkably, a similar dynamics is also observed for 
initial conditions with arbitrary initial angles between 
the angular momenta Li and L 2 . We start again from 
the linear 3D harmonics oriented this time arbitrary with 
respect to each other and instantaneously switch-on the 
nonlinear interaction. Examples of the vortex precession 
and nutation for different initial [3 in @ are shown in 
Fig. 4(a). The dynamics obtained for (3 < 7r/2 degrees 
are clearly different from those for /3 = 7r/2 degrees since 


FIG. 4: (Golor online) (a) Precession of nonlinear vortices 
grown from the linear harmonics with N = 50.11 and various 
values of [3. (b) Dependence of the precession period on the 
initial angle /3. 

for (3 'K/‘I instead of spikes smooth oscillations appear 
which correspond to a precession with a forward release 
speed EH- For small /3 the amplitude of the nutations 
decreases significantly and is hardly observable. Also for 
smaller values of (3 precession is much faster, therefore the 
opening angles of axes cones decrease with (3 as we can 
see comparing the scales of panels of Fig. 4(a) . Fig.|^b) 
shows the dependence of the full precession period on the 
initial angle f3 (for the sake of compactness we merged 
motion of two axes into single panels showing evolution of 
axes of ipi and ^2 with blue and red colors respectively). 
For all considered N the period decreases strongly when 
P deviates from 7r/2. 

To conclude, we have reported that a two-component 
BEG with repulsive inter- and intra-species interactions 
can support very robust non-coaxial vortices, with ap¬ 
proximately orthogonal vortex lines in each of the com¬ 
ponents. These cross vortices can be excited from nonin¬ 
teracting eigenstates of the 3D linear harmonic oscillator 
by sudden switch-on the nonlinearity and they are char¬ 
acterized by remarkably persistent dynamical regimes of 
precession with nutation resembling the motion of a rigid 
body. The precession is generated by the initial condi¬ 
tions, while the nutation occurs due to the force caused 
by the inter-species interactions and is associated with 
deformations of vortex shapes. 
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